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Abstract 

We study multiple elastic collisions of a block and a ball against a rigid wall in one dimension. 
The complete trajectory of the block is solved as an analytic function of time. Near the turning 
point of the block the force carried by the ball is proportional to 1/x 3 , where x is the distance 
between the wall and the block, in the limit that the block is sufficiently heavier than the ball. 
This is a unique pure mechanical system that reveals dipoledike repulsion. 



I. INTRODUCTION 



An elastic collision is a useful mechanical problem with which one can study energy, linear- 
momentum, and angular- momentum conservation in an explicit way. 1 While the final state 
is uniquely determined in a two-body elastic collision in one dimension, there are multiple 
solutions if three or more particles are involved. The reason is that the two constraints 
from energy and momentum conservation are not enough to determine the three or more 
final-state momenta. For example, a unique solution for a Newton's cradle of three or more 
pendula is obtained only if any two neighboring pendula are separated so that each collision 
involves only two pendula. 

Various studies on such chain collisions have been carried out for the instantaneous con- 
tact force, 2-7 for the spring force, 8-11 and for the contact Hertz force. 7 ' 8 ' 10 ' 12 In addition, there 
are numerous studies of the collisions of aligned balls in one dimension. 13-21 Harter investi- 
gated the velocity gain in such a chain collision. 16 Patricio studied the effect of the Hertz 
contact force in detail. 7 Redner considered the one-dimensional collisions of two cannonballs 
with a table-tennis ball sandwiched between them to derive a simple relation between that 
elastic collision and a corresponding billiard system. 22 

In this paper, we consider the multiple collisions of two particles in one dimension against 
a rigid wall. As is stated above, this problem is, in principle, exactly solvable and we attempt 
to find the analytic solution to this problem. Many people, who have played or watched 
table-tennis games, are familiar with this multiple collision. When the server is waiting for 
the defender, the server often bounces the ball against the table repeatedly with the racket. 
While the player presses down the ball with the racket, the speed of the ball keeps increasing 
and one can hear an interesting high-frequency sound near the turning point of the racket. 

Whelan et al. studied a similar case under gravity focusing on its chaotic behavior. 23 
Because the main feature of the problem stems from the velocity amplification of the light 
target against the heavy incident particle, we ignore any external forces other than the 
contact interaction between any two colliding objects. We construct a simple model system 
consisting of a block and a rigid wall that sandwich a ball. We derive the analytic expression 
for the complete trajectory of the block as a function of time. This new analytic solution 
reveals that, when the ball is sufficiently lighter than the block, the ball mediates a force 
proportional to 1/x 3 near the turning point of the block. Here, x is the distance between the 
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wall and the block. This emergent dipole interaction is unique in a pure mechanical system. 
In comparison with a previous study in Ref. 24 for a case with many balls, our derivation of 
the 1/x 3 rule is based on a complete analytic solution with a rigorous error analysis. 

This paper is organized as follows. In Sec. II we introduce our model, which simplifies the 
multiple collision problem of the table-tennis ball. In Sec. Ill we construct the simultaneous 
recurrence relation for the velocities of the colliding particles in each collision. Those for 
the position and time are also given. We determine the complete trajectory of the block as 
a function of time. In Sec. IV, we investigate the dynamics of the system in the continuum 
limit, where the block is sufficiently heavier than the ball and the block is near the turning 
point. Our conclusion is given in Sec. V, followed by the appendices, which contain useful 
mathematical formulas and some details of calculations. 

II. THE MODEL SYSTEM 

In this section, we define our model system. We restrict ourselves to elastic collision 
in one spatial dimension and ignore any external forces other than instantaneous contact 
forces. 

As shown in Fig. 1, the model system consists of a block with mass M, a ball with mass 
m, and a rigid wall at x < 0. The unit vector along the positive x axis is given by x and 
the mass ratio is defined by 



At time t — t\ — the block hits the ball at x — x\ — L, where initial velocities of the block 
and ball are 



respectively. After the first collision, the ball bounces against the wall and then they make 
the next collision. We call P n (t n ,x n ) the n-th collision point between the ball and the 
block, where t n and x n are the time and position. The velocities of the block and the ball 
right after the n-th collision are U n and u n , respectively. The total number of collisions N 
is the smallest n that satisfies the condition U n > \u n \. We also define the time interval 
At n = t n+ i — t n between P n and P n+ \. Useful relations involving the computation of t n 
using At n are given in Appendix A. 



a = m/M. 



(1) 



Uq = —U and u = 0, 



(2) 
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The elasticity of the collisions requires the conservation of energy as well as the linear 
momentum in the n-th collision: 



MU n - mu n = MU n+ i + mu n+ i, (3a) 
l -MU 2 n + l -mu 2 n = \MUl +l + l -mu 2 n+l . (3b) 

Note that the sign of u n is flipped in the first equation because the ball has bounced against 
the wall. With the initial conditions in Eq. (2), one can solve U n and u n recursively The 
first two pairs of the solutions are 





(4a) 

U I 1 - 6a + a 2 , 

(4b) 



(1 + a) 2 I 4(1 -a) 

Based on Eqs. (4), we classify the ranges of the mass ratio a that determine the value of N. 

• If a > 3, then N = 1 because U\ > —Ui. The block bounces back after the first 
collision and they never collide again. At a — 3, U\ — —U\ — \ U. 

• If 1 < a < 3, then = 2 because Ui < —u\ and U 2 > u 2 > 0. At a = 1, Ui = 0, 
Ui = —U, U 2 = U, and u 2 = 0. 

• If a < 1, then N > 2. Especially for an a < 5 - 2y/E, N > 3. At a = 5 - 2^, 
C/ 2 = — U2- 

In the remainder of this paper, we consider the case < a < 1 that involves multiple 
collisions. 

III. ANALYTIC SOLUTIONS 

In this section, we determine the point P n (t n ,x n ) of the n-th collision between the ball 
and block, and the velocities U n and u n immediately after the collision. 
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A. Computation of U n and u 



We can reduce the recurrence relations for the velocities in Eq. (3) into the form 

U n +i i , i « i 

(5) 

u n +i 

where the 2x2 matrix A is defined by 




1 — a —2a 



1 + a I 2 1-a 
Applying Eq. (5) recursively, we determine U n and u n : 



(6) 





(7) 



Although our initial conditions are given in Eq. (2), the relation (7) is valid for any initial 
values of Uq and uq. 

The computation of U n and u n is rather tedious and we summarize the calculation in 
Appendix B. The results are 




-U 



( cos n6 ^ 
1 

— = sin n9 



(8) 



where the parameter 9 is related to the mass ratio a as 

tan | = y/a. (9) 

Let us interpret the results in Eq. (8). Because the total kinetic energy is the sum of quadratic 
functions of U n and u n , the conservation of the total kinetic energy implies that U n and u n 
can be parametrized as a cosine and a sine of a common phase. The phase increases as a 
multiple of an elementary phase 6 because the two eigenvalues of the non-Hermitian matrix 
A are a unimodulus complex number e ld and its complex conjugate. According to these 
results, the block slows down to reach the turning point around n6 ~ |, where \u n \ becomes 
the maximum. The maximum value for \u n \ is larger than the initial speed U of the block 
with the enhancement factor l/y/a, if < a < 1. From that moment to the region n6 pa ir, 
the block bounces back and the ball loses the kinetic energy. If n — iV, then the velocity 



— u n of the ball after the bounce against the wall must satisfy the constraint U n > —u n . 
Therefore, N must be the smallest integer that satisfies 

U N + u N >0 and N > (10) 

B. Computation of x n 

The recurrence relation of x n is derived in Eq. (A2) of Appendix A: 

— = (u n -U n )/(u n + U n ). (11) 

x n 

By making use of Eqs. (8) and (9) and the elementary trigonometric identity sin a cos (3 ± 
cos a sin (3 = sin(a ± (3), we find that 

x n +i _ sin(n - \)6 



x n sin(n + \)9 

With the initial condition x\ = L, we find the analytic expression for x n as 



(12) 



n-l 
k 



n x k+i Lsin g 

i x k -sm(n-l)d- {l6) 



The analytic expression for x n is new. From this compact expression, it is easy to see that 
Xn decreases to reach the minimum value ~ |(| — 1), where the denominator 

reaches its maximum. The turning point x m i n is reached at n = Ni and the last collision is 
made at xn- 



C. Computation of N\ and N 

From Eq. (13), we can determine that Ni is 

N 1 = 



TT 



291 + ' < 14 > 

where [x} + is the smallest integer that is greater than or equal to x. For a small a, 9 2y / a 
and therefore N± In analogy to the derivation of Eq. (13), we can simplify U n + u n 

as 

U n + u n = \ e ■ ( 15 ) 

sin | 



6 



Substituting Eq. (15) into Eq. (10), we find that 



N 



7T 



2v^ 



(16) 



where the final expression is valid for small a. This result is consistent with Eq. (6) of Ref. 22 
with the uncertainty ±1. The quantity N max of that reference is approximately equal to 2N 
because, in that reference, the collision between the ball and the wall is also counted. 

In Fig. 2, we show the total number of collisions iV as a function of a. The results in the 
figure are consistent with classification given at the end of Sec. II. The discontinuities of N 
are at a = a.^ where ctj is the minimum value of a that satisfies N = i. The first two values 
of (Xi are ct\ = 3 and a 2 = 5 — 2\/5. 



D. Computation of t n 

We have determined U n , u n , and x n for all n — 1, • • • , N . In this section, we use these 
results to compute the collision times t n . The computation of the time interval At n is given 
in Appendix C. By making use of Eq. (A4), we sum the time intervals At n in Eq. (C3) to 
find 

tan(n — \)9 

where r is the time interval for the block to reach the wall when the ball is absent [a — 0): 

t = L/U. (18) 

This analytic expression (17) for t n is new. 
The complete trajectory of the block is 

x(t)=x n + U n (t-t n ), t n <t<t n+1 , (19) 

where x n , U n , and t n are given in Eqs. (13), (8), and (17), respectively. In Fig. 3, we plot the 
trajectory x(t) of the block as a collection of straight line segments. The vertices represent 
the collision points P n (t n ,x n ). The point on the vertical axis represents the initial collision 
point Pi(0, L). The turning point is P/v^mm, ^min)- When a = 0.25 (0.1), we have N± — 2 
(3), and N — 3 (5). The x-t plot is convex-downward because every impact on the the block 
is along the positive x axis. 



(17) 
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E. Computation of x m i n and i m i n 

Next we find the turning point x min — x^ and corresponding time t min — t^. According 
to Eq. (14), (Ni — \)0 is bounded by 

\(*-0) < {Ni-\)6< l(n + 9). (20) 

Therefore, cos | < sin(Aq — |)0 < 1. Substituting this constraint into Eq. (13) with n = Aq, 
we can determine the range of rr min as Lsin| < < Ltan|. With the values for sin | 
and tan \ in Table I, we obtain 



2 

L\fa 



< x mhl < L^. (21) 



v / T+ a 

When a is small, x min ps L^fa. From the range in Eq. (20) we find that cot | < tan(A r 1 — 
\)6 < oo or — oo < tan(iVi — \)9 < — cot |. Therefore, considering both cases, we obtain 

r(l-a) < t mi „ < r(l + a), (22) 

where we used tan 2 \d = a. The errors of x m i n and t min are ±.\L^/a{l — 1/y/l + a) and 
±ar, respectively. 

It is interesting to notice that when a is sufficiently small, the time interval i m i n for the 
block to reach the turning point is approximately the same as r, which is the time spent for 
the block to hit the wall when the ball is missing. 



IV. CONTINUUM LIMIT 



According to Eq. (CI), At n becomes small when a is sufficiently small and n ~ A r 1 , where 
sin(n — |)0 ~ 1. This does not require that a — > + . We call this the continuum limit. In 
this limit, both x n and t n can be treated to be continuous variables. 

In this section, we investigate the dynamics of the system in the continuum limit. As the 
first step, we express the collision point x n in terms of t n . Next, we introduce a differentiable 
function x(t) that may represent the trajectory of the block in the limit. By making use of 
the kinetic energy conservation, we interpret the kinetic energy of the ball as the potential 
energy of the block. Expressing the kinetic energy of the ball in terms of the differentiable 
function x(t), we compute the force on the block in the continuum limit. 
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A. Relation between x n and t 



We observe that both x n in Eq. (13) and t n in Eq. (17) depend on n explicitly. In this 
section, we find the analytic expression for function of t n by eliminating the explicit 

n dependence. 

From Eqs. (13) and (17), we have 

x n cos(n— \)6 = L ■ - ^^ Zr ? (23a) 

Vl + a 

x n sm(n-\)6 = L \J y~^' ( 23b ) 

where sin | and tan | are replaced with the values listed in Table I. By making use of the 
fact that sin(n — \)6 > for all n < N, we can express x n as a function of t n : 

L 



x r 



^{l-tjrf + a, (24) 



where we used relations in Eq. (23) and Table I. 

In Appendix D, we have shown that the complete trajectory x(t) of the block is always 
bounded by 

fit) < x(t) < g(t), (25) 
where the different iable functions f(t) and g(t) are defined by 

fit) = -±=y/(l-t/T)i + a, (26a) 

V 1 + OL 

g(t) = L^il-t/rf + a. (26b) 
As shown in Fig. 3, the lower bound f(t) passes every collision point P n and the upper 



bound g(t) is tangent to every line segment P n P n+ \. 

In the continuum limit, it is convenient to use the arithmetic average x(t) = \[g{t) + f(t)] 
to represent the trajectory of the block with the uncertainty Sx(t) = \[g{t) — f{t)\. 

x(t) = ML^il -t/r) 2 + a, (27a) 

5x(t) = ALy/(l-t/T) 2 + a, (27b) 

i i 
where M = \[l + (1 + a) 2] and A = \[l - (1 + a) 2]. In the limit a -> + , all of the 

curves x(t), f(t), and git) collapse into x(t). In addition, 

lim x(t) = L\l-t/r\, (28) 



which is the case that the ball is absent. 

In the continuum limit we can compute the velocity and the acceleration of the block 
using x(t): 

MUil-t/r) 



x(t) 
x(t) 



(29a) 
(29b) 



[(1 -t/r) 2 + a] 3 / 2 ' 

The values for x, x, and x must be understood as the time-averaged values. The uncertainties 
of x(t) and x(t) can be obtained by replacing M with A. 



B. Kinetic energy of the block 

Based on Eqs. (29a) and (27a), we can determine the kinetic energy K of the block in 
terms of t or x as 

M 2 {l-t/r) 2 



Kit) = K - 



K(x) = K J\f- 



l-t/r) 2 + a 
aU 2 L 2 



1 - 



(30a) 
(30b) 



where K = \MU 2 is the initial kinetic energy, which can be approximated as K K(t = 
— oo). For a small a, the relative error of the kinetic energy is ±|ct. 

C. Potential energy of the block 

We recall that the total kinetic energy is conserved during the whole process. Therefore, 
the sum of the kinetic energies of the block (K) and the ball (<&) must be 



K + $ = K . 



(31) 



We interpret the kinetic energy $ of the ball as the potential energy of the block. By making 
use of Eqs. (30a), (30b), and (31), we find that 

M 2 {l-t/T) 2 



$(x) 



K 
K 



1 - 



(1-t/r) 



a 



1-M 2 + 



aM A L 2 



(32a) 
(32b) 
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Next we compute the force, 



F(*) = -^ = «^.£ (33) 



The force on the block has the same position dependence as an electric dipole interaction: 
F(x) ~ 1/x 3 . In addition, x = x min when t = r. Therefore, the duration of time to reach 
rr min is the same as the case when the ball is absent. As is stated earlier, the derivation of 
the 1/x 3 rule originates from the analytic solution in Eq. (27), which is in contrast to the 
approach in Ref. 24 . 

We can compute £ mm in the continuum limit. Setting $(a; mm ) = K , we find that 

•^min 

= MLy/a, (34a) 
t min = r, (34b) 

These values are consistent with Eqs. (21) and (22) within errors. 



V. CONCLUSION 



We have considered the one-dimensional elastic collisions of a ball and a block against a 
rigid wall. The initial state of the ball is u = at x = L and that of the block is U = —U. 
The trajectory x(t) of the block is completely determined in an analytic form. The analytic 
expression for the total number N of collisions between the block and the ball was also 
derived. The turning point of the block is ~ L^fa at t ~ r = L/U, which is the time 
interval for the block to hit the wall when the ball is absent. Here a is the ratio of the mass 
of the ball to the mass of the block. 

In the continuum limit where a is small and x(t) is near the turning point, the trajectory 
can be approximated as a differentiable function x(t) in Eq. (27a). Because the total kinetic 
energy is conserved in this system, one can think of the kinetic energy of the ball as the 
potential energy of the block. Based on this idea, we have computed the force on the block 
and found that the force is proportional to 1/x 3 . It is remarkable that this is a unique pure 
mechanical system that reveals repulsive dipole-like interaction. 
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Appendix A: Useful Formulas to compute t n 

The time interval between P n and P n+ i is defined by At n = t n +i — t n which can be 
expressed in terms of U n or u n as 

At n = Xn+ \~ Xn = ^L, (Al) 

where s n = x n + x n+ i is the distance that the ball travels between P n and P n +i- Assuming 
u n < 0, we find that 

x n+ i = x n (u n - U n )/(u n + U n ). (A2) 
Substituting Eq. (11) into Eq. (Al), we can express At n as 

At n = --^. (A3) 

The n-th collision time t n can be computed as 

71-1 



t n = ^At fe . (A4) 



k=i 



Appendix B: Computation of U n and u r 

In this section, we solve Eq. (7), 



VA -AT "'], (Bl) 



M \ WO 

where the matrix A is defined in Eq. (6), 



1 + a \ 2 1 - a 
We make a transformation diagonalizing A such that 



1 / 1 - a -2a , 
A = TTT\ • ( fi 2) 









| , R^AR 


I J«n / 




I J«n / 





Ai 
A 2 



(B3) 



where i? is a 2 x 2 matrix and Aj's are the eigenvalues of the matrix A. Then, U n and u n are 

V„\ = AlWo\ [ W °\ = R' 1 [~ U \ . (B4) 
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The i-th column vector i?; of R satisfies the equation, 

AR t = X t Ri, (B5) 

for i — 1 and 2. The corresponding secular equation det[A — Ajl] = 0, where 1 is the 2x2 
identity matrix, is 

a(l + A,) 2 + (l-A,) 2 = 0. (B6) 

Solving this equation, we obtain 

X, = e +w , R 1= I " I , (B7a) 



-i/y/a 

a, = i -r /;. ( 1 | . (Bib) 

where the parameter 9 is related to the mass ratio a as 

tan | = y/a. (B8) 

In Sec. IV of Ref. 23 , the authors considered the collisions of two balls near the ground, which 
can be compared to our results by neglecting gravity. The matrix A in Eq. (6) is equivalent 
to the matrix M in Eq. (13) of Ref. 23 . The elements in the second column have opposite 
signs because of the difference in the definition of u n . The eigenvalues \ in Eq. (B7) and 
the parameter 6 in Eq. (9) are consistent with those in Eq. (15) of Ref. 23 . 

Here, we consider only the case < a < 1. Note that 9 = for a = 0, < 9 < | for 
< a < 1, and 9 — | for a — 1. In the second and third columns of Table I, we list the 
values for trigonometric functions for \9 and 9, respectively. Now we determine R and R~ l 
as 

r= i 1 , ,\ I i- - .; i ^ i (bo) 




-i/y/a i/\fa J \ 1 

Substituting Eq. (B9) into Eq. (B4), we finally determine the velocities U n and u n as func- 
tions of n and a: 



Mr. 



cos n9 

1 . J- (BIO) 
sin n9 

a 
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Appendix C: Computation of the time interval At 



The time interval At n = t n+1 — t n can be computed by substituting x n in Eq. (13) and 
U n + u n in Eq. (15) into Eq. (A3). 



At 



2t sin 



" " sin(n-|)0sin(ri+§)0' 



(CI) 



where r is the time interval for the block to reach the wall when the ball is absent (a — 0), 
that is, r = L/U. Substituting a = (n—\)0 and (5 = (n+\)9 into the following trigonometric 
identity, 



sin(/3 — a) 



sin a sin , 



we can simplify Eq. (CI) as 

At n = r 



tan 



= cot a — cot /3, 



tan 



tan(n — -AO tan(n + AO 



(C2) 



(C3) 



Appendix D: Determination of the functions f(t) and g(t) 

The complete trajectory of the block x(t) is a set of line segments that connects consec- 
utive collision point such that x(t n ) = x n , where 

L 



X j 



VI + a 

We introduce a differentiable function, 

fit) L 



^(l- t n /rf + a. 



(Dl) 



V(l-*/r) 2 + a, 



(D2) 



that passes every collision point P n and is convex downward. Therefore, at any time t G 
[0, £jv], /(t) < x(t), where the equality holds only at t — t n for n — 1, • • • , iV. 
The time derivative of Eq. (D2) is 

U 1-t/r 



/(*) = 



(D3) 



+ a 

By making use of the special values of the trigonometric functions in Table I, we find that 

/(*„) = -U cos §cos(n-|)0 = ±(U n -! + U n ), (D4) 
where we have used cos(a + /3) + cos(o; — 0) — 2 cos a cos /3 for a = (n — |)0 and /3 = 

14 



If there exists a curve g(t) that is tangent to every line segment P n P n+ \ at t' n G [t n ,t n+ i] 
and convex downward, then g(t) is an upper bound of x(t) for all t G [0, tjv]- We require 

9(t'n) = U n = -U cos n9, (D5) 

and we set x' n = g(t' n ). According to Eq. (D4), g(t' n ) = f(t n+ i)/ cos |. Therefore, we find 
that 

U(l - t/r) _ . 

9(t) = ~ , K , (D6a) 

a/(1 - t/r) 2 + a 

t' n = t i =T(l--v/acotra0), (D6b) 

n+ 2 

. . . L sin f , 

< = = ita^f ■ (D6c) 

Note that t n < t' n < t n+1 and g(t n ) = —U cos (n — |) 0. Integrating in Eq. (D6a) over t 
and imposing the boundary condition (D6c), we obtain 



g(t) = Ly/(l-t/ry + a. (D7) 

As a result, the trajectory x(t) is bounded by 

f(t)<x(t)<g{t), (D8) 

for any t G [0, £jv]. The left and right equalities hold at t — t n and at t — t' n = t n+ i, 
respectively, for all n. In Fig. 3, we show the trajectories of the block x(t) at a = 0.25 and 
a — 0.1. As shown in this figure, the lower bound f(t) passes every collision point P n and 



the upper bound g(t) is tangent to every line segment P n P n+ i. 
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Tables 



f(x)\x 



cosx 



sin a; 



tanx 



fa 



1 — a 
1 + a 

2Ja~ 



1 - t n /T 



V(l-t n /r) 2 + a 
fa 



1 + a 
2 



V(l-Wr) 2 + a 



a 



o 

1 — a 



a 



1 - tn/r 



TABLE I: The values for cosx, sinx, and tanx for x = 5$; Q, and (n — \)0. The values at x = n6 
and (n + can be obtained by using the identities such as cos(x ± 5) = cos x cos 5 =F sin x sin (5 
and sin(x ± 5) = sin x cos 5 ± cos x sin 5. 
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FIG. 1: The initial condition of the model system. 
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FIG. 2: The total number of collision N as a function of the mass ratio a. qi = 3 and CK2 = 5 — 2y5 
are the minimum values of a to have N = 1 and N = 2, respectively. Asa->0 + ,iV-> ti^Ta — 2] + 
which diverges to infinity. 
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FIG. 3: (a) The collision points P n (t n ,x n ) at a = 0.25, where N\ = 2 and N = 3. The trajectory 
of the block x(t) consists of line segments that connect every collision point P n . The lower bound 
f{t) passes every collision point P n and the upper bound g(t) is tangent to every line segment 
P n P n +i. The continuum representation x{t) of the trajectory is the arithmetic average of git) and 
fit), (b) The same as (a) except that a = 0.1, N\ = 3, and N = 5. 
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